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ABSTRACT 

We report the results of a search for pure rotational molecular hydrogen emission from the cir- 
cumstellar environments of young stellar objects with disks using the Texas Echelon Cross Echelle 
Spectrograph (TEXES) on the NASA Infrared Telescope Facility and the Gemini North Observatory. 
We searched for mid-infrared II2 emission in the S(l), S(2), and S(4) transitions. Keck/NIRSPEC 
observations of the II2 S(9) transition were included for some sources as an additional constraint on 
the gas temperature. We detected H2 emission from 6 of 29 sources observed: AB Aur, DoAr 21, 
Elias 29, GSS 30 IRS 1, GV Tau N, and HL Tau. Four of the six targets with detected emission are 
class I sources that show evidence for surrounding material in an envelope in addition to a circum- 
stellar disk. In these cases, we show that accretion shock heating is a plausible excitation mechanism. 
The detected emission lines are narrow ('^10 km s~^), centered at the stellar velocity, and spatially 
unresolved at scales of 0.4", which is consistent with origin from a disk at radii 10-50 AU from the 
star. In cases where we detect multiple emission lines, we derive temperatures > 500 K from '^1 
of gas. Our upper limits for the non-detections place upper limits on the amount of H2 gas with T 
>500 K of less than a few Earth masses. Such warm gas temperatures are significantly higher than 
the equilibrium dust temperatures at these radii, suggesting that the gas is decoupled from the dust 
in the regions we are studying and that processes such as UV, X-ray, and accretion heating may be 
important. 

Subject headings: circumstellar matter, infraredistars, stars:planetary systems, protoplanetary disks, 
stars:individual(AB Aur, DoAr 21, Ehas 29, GSS 30 IRS 1, GV Tau N, HL Tau), 
stars:prc-main sequence 



1. INTRODUCTION 

Studying the structure and evolution of circumstellar 
disks is crucial to developing an understanding of the 
process of planet formation. Observations of dust emis- 
sion and modeling of the spectral energy distributions 
(SED) of disks have revealed much about the dust com- 
ponent frorn a few stellar radii out to hundreds of AU 
(|Zuckerman ll200iT). While circumstellar disks are com- 
posed of both dust and gas, the gas component dominates 
the mass of the disk, with molecular hydrogen (H2) be- 
ing the most abundant constituent. In order to develop 
a complete picture of the structure and evolution of pro- 
toplanetary disks, it is important to observe the gas. 

Observations at different wavelengths probe different 
disk radii. Submilli meter observations sa mple gas at 
large radii (> 50 AU) (ISemenoy et al.ll2005D. while near- 
infrared CO (jNajita et al.ll2003l : iBlake fc Boogertll2004( l 
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and H2O (|Carr et al.ll2004HThi fc ^ '20051 observations 
allow for study of the inner few AU. Spectral lines in the 
mid-infrared (5-25 /^m) provide a means to investigate 
gas in the g iant planet region of the disk and beyond 
(10-50 AU) (jNajita et al.ll2007a^ . 

Several mid-infrared spectral diagnostics have been 
shown to be useful p robes of gas in disks. These include 
[Nell ] at 12.8 ^m (jPascucci et all l2007t iLahuis et aP 
2007; Hcrczeg et al.' '2007|i^ H2O rotational transitions 
Carr fc Najito 2008; Salvk et all 12008 1. [Fel] at 24 /im 



Lahuis et al.ll2007l ). and, based on a theoretical analy- 



sis of debris disks, [SI] at 25.2 /im, and [Fell] at 26 /im 
(|Gorti fc Hollenbachl[200l . 

Molecular hydrogen should make up the bulk of the 
mass in disks, but is a challenge to detect. Bright far- 
ultraviolet (FUV) H2 emission from classical T Tauri 
stars may be produced in the irradi ated disk surface 
(|Herczeg et al.ll2002l : lBCT"gin et al.ll2004D . At Ion ger wave- 
lengths, rovibrational and pure rotational transitions are 
generally weak because H2 lacks a permanent dipole mo- 
ment. Near-infrared emission in the w=l-0 S(l) rovibra- 
tional transition of H2 has been detected from T Tauri 
stars (iBarv et al . 2003; Ramsav Howat fc Greavesll2007l : 
ICarmona et al.l 120081)1) and may be the result of excita- 
tion by UV and X-ray irr adiation (Nomura ct al. 200'3; 
iGorti fc HoUenbachl l2008h . Near-infrared adaptive op- 
tics fed, integral field spectroscopy of six classical T Tauri 
stars that drive powerful outflows has revealed that most 
of the H2 emission is spatially extended from the con- 
tinuum (Beck et al. 2008). The properties of the emis- 
sion are consistent with shock excitation from outflows 
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or winds rather than UV or X-ray excitation from the 
central star. The FUV H2 emission probes gas between 
2000 and 3000 K ()Herczeg et al.ll2004l: iNomura fc Millal 
I2005D and the 2.1 2 fim H2 line traces gas at T>1000 K 
( Bary et al.ll2003l ). The mid-infrared H2 lines considered 
in this paper are most sensitive to gas at lower temper- 
atures. 

Owing to their small Einstein A- values, the pure ro- 
tational mid-infrared H2 lines remain optically thin to 
large column densities (N//^ > 10^'^ cm~^) and will be 
in LTE at the densities found in disks. For a disk with 
a strong mid-infr'ared continuum, the dust becomes op- 
tically thick well before the H2 lines. Observable line 
emission is present only when there is a temperature in- 
version in the atmosphere of the disk or if there is a 
layer of dust-depleted gas separate from the optically 
thick dust. The process of dust coagulating into larger 
grains or settling out of the disk atmosphere can allow a 
larger column of gas to be observed. A disk that has an 
optically thin mid-infrared continuum, implying a very 
small amount of dust in the disk or dust grains that 
have grown large compared to mid-infrared wavelengths, 
would allow the entire disk to be observable. However, 
it is not known whether such disks have large quantities 
of gas and, in the absence of gas heating through col- 
lisions with dust grains, another heating mechanism is 
necessary in dust-free environments, such as UV or X- 
ray heating (iGlass gold et al.l [2OOI iGorti fc HollenbachI 
[2OOI INomura et al.ll2007ll 

A number of groups have searched for H2 emission 
from protoplanetary disks in recent years. iThi et al\ 
(|2001h reported the detection of several Jupiter masses 
of warm gas in a sample of disk sources based on In- 
frared Space Observatory (ISO) observations of the H2 
S(0) and S(l) lines. However, follow-up observations 
from the ground with i mproved spatia l resolution did not 
confi r m these results (iRichter et all 2002; S heret et all 
120031 ISako et all l2005f) . IRichter et al.. (2002) used the 
Texas Echelon Cross Echelle Spectrograph (TEXES) on 
the NASA 3m Infrared Telescope Facility (IRTF) to set 
upper limits on t he warm gas mass w ithin the disks of 
six young stars. iSheret et all (|2003f ) searched for H2 
emission using MICHELLE on the United Kingdom In- 
frared Telescope and set upper limits on the emission 
from disks around two stars. The first group to use 
an 8-meter cl ass telescope in t he search for molecular 
hydrogen was ISako et all ((20051) using the Cooled Mid- 
Infrared Camera and Spectrometer on the 8.2 m Sub- 
aru telescope to set upper limits for emission in the 
S(l) line around four young stars. Using the Infrared 
Spectrograph ( IRS) aboard the Spitzer Space Telescope, 
iPascucci et al.l (2006) reported the non-detection of H2 
lines in their samp le of 15 young Sun-like stars, while 
iLahuis et all (|2007[ ) detected the S(2) and S(3) Hnes in 
~8% of the 76 circumstellar disks in their sample. Re- 
cently, ground-based observations with high resolution 
spectrometers on 8-meter class telescopes have produced 
both detections of the mid-infrared H2 lines in the Her- 
big Ae stars AB Aur an d HD 97048 (|Bitner et all 120071 : 
iMartin-Zai'di et alll2007f) and stringent upper limits in a 
sample of six Herbig A e/Be stars and one T Tauri star 
(jCarmona et al.ll2008aD . 

Three mid-infrared pure rotational H2 lines are observ- 
able from the ground: S(l) (A = 17.035 fim), S(2) (A = 



12.279 fim), and S(4) (A = 8.025 iim). When multiple 
optically thin lines are observed, line ratios permit the 
determination of the temperature and mass of the emit- 
ting gas. Ratios of these three lines are most sensitive 
to temperatures of 200-800 K. Two additional pure ro- 
tational H2 lines are observable near 5 fim: S(8) (A — 
5.053 ^m) and S(9) (A = 4.695 liui), extending our tem- 
perature sensitivity to hotter gas. The high spectral res- 
olution possible with an instrument like TEXES (Lacy 
et al. 2002) increases our sensitivity to narrow line emis- 
sion and helps determine the location of the emission. By 
making observations at high spectral resolution, we max- 
imize the line to continuum contrast while minimizing 
atmospheric effects by separating the lines from nearby 
telluric features. A further benefit of high spectral reso- 
lution is that we are able to estimate the location of the 
emitting gas if coming from a disk under the assumption 
of Keplerian rotation. 

In this paper, we present the results of a search 
for molecular hydrogen emission in disk sources using 
TEXES on both the IRTF and Gemini North telescopes. 
We observed 29 sources spanning a range of mass, age, 
and accretion rate in order to constrain the amount of 
warm gas in the circumstellar disks of these stars. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed 29 sources using TEXES on both the 
NASA IRTF telescope during 2002-2005 and on Gemini 
North in 2006 and 2007 under program IDs GN-2006A- 
DS-3, GN-2006B-Q-42, and GN-2007B-C-9. The spec- 
tral resolution of the observations was >60,000 for the 
S(l) line and >80,000 for the S(2) and S(4) lines. Due to 
telluric absorption from water vapor close to the S(l) and 
S(2) lines, we observed these settings both when telluric 
water vapor levels were low and the Earth's motion gave 
an additional redshift towards the source. Removal of 
background sky emission was achieved by nodding the 
source along the slit and subtracting nod pairs. On the 
IRTF, the TEXES sht widths were 2" at the S(l) setting 
and 1.4" at S(2) and S(4). For our Gemini observations, 
the slit widths were 0.81" at S(l) and 0.54" at S(2) and 
S(4). 

We observed telluric standards at each setting for use 
as divisors to correct the spectra for atmospheric absorp- 
tion. Asteroids work well as telluric calibrators at 12 
and 17 /im, while we used early- type stars at 8 /im. We 
observed flux standards for some sources. Where possi- 
ble, however, we normalized the continuum level to agree 
with photometric measurements obtained from Spitzer, 
ISO, or ground-based observations. Tables [2l [H and [5] 
list references for the continuum fluxes. Data reduc- 
tion was carried o ut using the standard TEXES pipeline 
(|Lacv et al.ll2002D . 

Our sample was chosen to include sources with a range 
in age, accretion rate, and mass. Included among our 
targets are class I sources with both a remnant envelope 
and a disk, classical T Tauri stars with optically thick 
disks, T Tauri stars with optically thin disks or inner 
holes, stars with high accretion rates such as FU Ori and 
Z CMa, Herbig Ae stars with more massive disks, as well 
as the debris disk around the star 49 Ceti. The combina- 
tion of TEXES with Gemini has allowed us to extend our 
survey to sources with mid-infrared continuum levels of a 
few tenths of a Jansky. Since the mid-infrared continuum 
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comes from warm dust grains, this aUows the inclusion of 
sources without much dust where larger columns of gas 
may be observable as long as there are additional heating 
mechanisms beyond gas/dust heating. Properties of the 
targets in our sample arc listed in Table [T] 

Observations of the H2 S(9) line for a subset of our 
sample were carried out between 2000 and 2005 using 
NIRSPEC on the Keck telescope as part of an ongo- 
ing M-band survey (iBlake fc Boogertll2004D at a spectral 
resolution of R=:25,000. The S(9) observations were car- 
ried out using a slit width of 0.43". For a description 
of the data reduction process for these observations see 
iBoogert et all (|2002bl ). 

3. RESULTS 

Table [2] summarizes our line detections. We detected 
H2 emission from 6 of 29 sources observed: AB Aur, 
DoAr 21, Ehas 29, GSS 30 IRS 1, GV Tau N, and HL 
Tau. Fluxes, line widths, and centroids were determined 
by fitting each line individually with a Gaussian pro- 
file, which describes the detected lines reasonably well in 
most cases. In four of the five sources observed from both 
the IRTF and Gemini, the line fluxes or upper limits from 
Gemini are smaller than those measured from the IRTF. 
Since the TEXES sht is larger on the sky at the IRTF, 
this suggests the H2 emission may be spatially extended. 
Figures [3]l9] show that the lines are all centered near the 
stellar velocity suggesting the II2 is associated with the 
targets. In all cases, the mid-infrared lines are narrow 
with observed FWHM near 10 km s~^. Correction for 
the TEXES instrumental line width of 4-5 km s~^ would 
give source line widths '-^l km s^^ smaller than observed. 
For the Keck/NIRSPEC S(9) observations listed in Ta- 
ble [U removal of the instrumental line broadening would 
give line widths 3-4 km narrower than observed. Al- 
though our spectra cover ±100 km s^^ from the stellar 
velocity for the S(2) line, and a similar amount to the 
red of the S(l) line, we see no evidence for emission at 
large Doppler shifts. We most often obtain the high- 
est signal-to-noise detections of the S(2) line, where the 
atmospheric transmission is higher and TEXES is more 
sensitive. 

For emission from optically thin gas, the line flux is 
given by 

NuX Aui xhv 

where is the number of II2 molecules in the upper 
state, V is the line frequency, and d is the distance to the 
source. The gas mass and temperature enter the equa- 
tion determining the number of molecules in the upper 
state, which, under the assumption of LTE at a single 
temperature, takes the form 



HLTau 
T . 465 K 



M ^ , , hcB 

X gN{2Ju + Ij^TTT^e 



2kT 



(2) 



where M is the H2 gas mass, m//^ is the mass of an II2 
molecule, is the nuclear statistical weight (3 or 1 for 
ortho or para) , B is t he rot ational constant for II2 taken 
from iJennings et al.l (Il984f). and E„ is the upper state 
energy taken from iMandv fc Martinf (|1993[ ). 

In Table [3] we show the derived temperature and mass 
for the emitting gas in each source based on both sin- 
gle temperature and two-temperature LTE models. Fig- 
ure [T] shows excitation diagrams for the three sources in 




1000 2000 3000 

Upper State Energy (K) 

Fig. 1. — Excitation diagrams for the three sources in our sample 
where we have observations of all three mid-infrared H2 transitions 
and detections of at least two. The points are based on Gaussian 
fits to each of the transitions and are plotted with l-cr error bars. 
The overplotted lines show the best-fit single temperature. 



our sample where we have observations of all three mid- 
infrared II2 transitions and detections of at least two. We 
assumed optically thin LTE II2 emission and constructed 
synthetic spectra for a range of temperatures and masses 
using the Gaussian line parameters determined from fits 
to each line individually. We simultaneously fit the syn- 
thetic spectra to all of the data for each source, including 
non-detections, and determined the best fit by minimiz- 
ing the square of the residuals. The errors listed in Ta- 
ble [3] are la based on the contour plot of the values. 
We present the results of our best fitting model in the 
appropriate figures, but defer discussion of the individual 
sources until the next section. 

Tables 21 and list our derived 3(t line flux upper limits 
for our IRTF and Gemini observations in the cases where 
no line was detected. The standard deviation {a) of the 
line fluxes was computed by looking at the distribution 
of values found when assuming a FWHM comparable 
to our line detections (10 km s~^) and summing over 
the number of pixels corresponding to that FWHM in 
regions of the spectrum with comparable atmospheric 
transmission. 

In the case of TW Hya, we used a FWHM of 5.5 km s^^ 
to compute the line flux upper limits based on the mea- 
sured width of the (2cr) S(2) feature. Figure [2] shows ob- 
servations of TW Hya made during the TEXES/Gemini 
engineering run in February 2006. There is a hint of a 
feature near the S(2) position but not a clear detection. 
Figures [131 pj [HI and [H show 3(7 S(l) and S(2) line 
flux upper limits assuming a FWHM for the lines of 10 
km s"^ overplotted on the continuum at the stellar ve- 
locity. The gas mass upper limits listed in Tables [4] and 
[5] were computed under the assumption that the H2 is 
in local thermodynamic equilibrium (LTE) and that the 
emission is optically thin. 

At the temperatures to which we are sensitive, T > 
200 K, the LTE value of the ortho/para ratio (OPR) 
is 3, which we assume for our calculations. However, 
nonequilibrium values of the OPR have been observed 
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AB Aur S(9) - 
Kcck/NIRSPEC - 



\Vi,o [km s-'l 

Fig. 2.— TW Hya TEXES/Gemini S{2) and S(l) overplotted 
with a Gaussian fit to the {2a) bump near the S(2) position and 
the 3(7 upper limit at S(l). The measured FWHM of 5.5 km s~^ 
for the Gaussian fit at S{2) was assumed to obtain the line flux 
upper limit at S(l). The S{1) and S(2) upper limits are listed in 
Tabl e [5] The dashed lin e indicates the systemic velocity of TW 
Hya l IKastner et al.|[T999D . 



in gas at these temperatures. Using ISO, iNeufeld et al.1 
(|1998f) derived an OPR of 1.2 towards HH 54 from gas at 
T ~ 650 K. They argued that the warm, shocked H2 gas 
they observed acquired its OPR at T < 90 K and that it 
has not had time to equihbrate to the L TE value at the 
higher temperature. iFuente et all (|1999D found an OPR 
between 1.5 and 2 from 300-700 K gas in a photodissoci- 
ation region (PDR) , which is of interest since the surface 
layers of the disks we may be obse rv ing have simi l arities 
with PDRs (Jonkheid et al.ll2004f) . iBitner et al.l (|2007f ) 
claimed that the surface brightness of the H2 emission 
observed towards the disk source AB Aur was similar 
to the Orion bar PDR, lending support to the idea that 
PDRs and disk surfaces have similar qualities. This was 
the result of a computation error. In fact, the surface 
brightness of the H2 emission from the Orion bar PDR 
is significantly larger than that of the AB Aur emission. 
We note that FUV pumping can lead to an apparent 
OPR in this range even in gas with an equilibrium OPR 
of 3 due to ortho-H2 pumping rates being reduced by 
self-shielding (jSternberg fc Neufeldl I1999D . If the OPR 
is actually less than 3 in the sources in our sample, a 
single-temperature fit to our observations would lead us 
to derive a higher gas temperature because the lowest 
energy transition we observe is the ortho S(l). 

As expected, the upper limits are generally more strin- 
gent for the Gemini data than for the IRTF data. The 
derived 3a line flux upper limits for our Gemini data are 
in the range of 10~^^ to 10"^"* ergs s~^ cm~^. Observa- 
tions with IRS on Spitzer of sources with low continuum 
fluxes give 3a line flux upper limits of 10~^^ to 10~^ ^ 
ergs s-i cm-2 (jPascucci et al.ll2006l : iLahuis et al.ll2007t ). 
The very low upper limits possible with Spitzer are for 
sources with continuum fluxes too low to be detected 
from the ground at the high spectral resolution available 
with TEXES. In one case where we have a source in com- 
mon, 49 Get, our derived upper li mits are nearly id entical 
to those obtained with Spitzer (jGhen et al.ll2006[) . Our 
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Fig. 3.— AB Aur NIRSPEC/Keck S(9) and TEXES/Gemini 
S(4), S(2), S(l) data from November 2006 observations overplotted 
with 2-component model fit. The dashed line shows the stellar 
velocity (Thi et al. 2001]). The increased noise in the S(l) spectrum 
blueward of the position of the S(l) line is caused by a telluric 
feature. 

upper limits constrain the amount of warm gas with opti- 
cally thin dust to be below several tens of Earth masses 
for a temperature of 200 K and less than a few Earth 
masses at temperatures above 500 K. 

4. DISCUSSION 
4.1. Individual Sources with H2 Detections 
4.1.1. AB Aur 

AB Aur is a Herbig Ae star surrounded by circum- 
stellar material extending to at least r ~ 450 AU 
(|Mannings fc Sargent ! 119971). It has a spectral type 
AO (iHernandez et a l. 2004) and is located 144 pc away 
(|van den Ancker et al. 1998 ). The disk surround ing the 
2.4 M© central star (van de n Ancker et al.]|1998f ) has an 
estimated mass of 0.013 Mq and an inclination of 17^3 
deg as determined by mode ling molecular line em ission 
at millimeter wavelengths (ISemenov "eran[200l . Ob- 
servations bv lGhen fc Jural ( 20031 ) show that AB Aur is 
variable i n the mid-in frared. 

iBitner et al.l ()2007[ ) described observations of the 
S(l), S(2), and S(4) H2 lines using TEXES. A single- 
temperature LTE model flt to the lines yielded a temper- 
ature of T = 670 K and M = 0.52 M© for the emitting 
gas. In this paper, with the addition of observations of 
the S(9) line, gas at a single temperature no longer fits 
the data. The two-temperature fit shows that essentially 
all the S(9) flux comes from a very small amount of gas 
(« 0.075 M0) at a temperature likely close to the dust 
sublimation temperature (Figure [3]). It also appears that 
the high temperature component is slightly blue shifted 
relative to the gas responsible for the lower energy lines. 
We observed AB Aur from Gemini on two separate occa- 
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Fig. 4.— AB Aur NIRSPEC/Keck S(9) and TEXES/Gomini 
S(4), S(2), S(l) data from October 2007 observations overplotted 
with the 2-component model fit derived using November 2006 data. 
The S(l) and S(2) lines are consistent with our November 2006 
observations, however the S( 4) line appears weaker. The dashed 
line shows the stellar velocity IIThi et al.l20011) . The increased noise 
in the S(l) spectrum blueward of the position of the S(l) line is 
caused by a telluric feature. 

sions, once during 2006 November and again in 2007 Oc- 
tober (Figures [3] and 2]) . In Figure 3] we show the obser- 
vations from 2007 overplotted with the two-temperature 
model based on the 2006 data. While the model is a 
reasonable fit to the S(I) and S(2) Hues, the S(4) line ap- 
peared weaker in our 2007 observations. Table [2] shows 
details of individual Gaussian fits to each of the lines from 
both years. The S(2) line fluxes are in excellent agree- 
ment while the S(I) and S(4) observations show differ- 
ences. It is possible that the source varied between the 
two observations. However, the agreement among the 
two different S(2) observations combined with the fact 
that we are more sensitive there than at S(f ) or S(4) casts 
doubt on this possibility. 

Modehng of the FWHM of the S(2) line profile account- 
ing for Keplerian, instrumental, and thermal broaden- 
ing suggests that the emission arises near 18 AU in the 
disk (Bitner et al. 2007). Spatially resolved mid- infrared 
continuum ir nages at 11.6 urn usi ng Michelle on Gem- 
ini North by iMarihas et al.l (|2006f) reveal a source size 
of 17 AU for the emission that remains after the sub- 
traction of a comparison PSF st ar. Comparis on with the 
passive flared disk model of DuU emond et al.l (.2001) sug- 
gests that, as the disk emerges from the shadow of the 
inner rim near 10 AU, dust grains in the surface layer 
are heated by direct stellar radiation, which produces the 
mid-IR continuum flux. Their derived average dust tem- 
perature in this region of 200 K is significantly lower 
than the gas temperature based on the II2 observations, 
implying an additional heating source for the gas. Likely 
candidates for the source of additional heating are X- 
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Fig. 5.— DoAr 21 TEXES/IRTF S(2) overplotted with Gaussian 
fit. The das hed line indicates t he systemic velocity of the p Oph 
cloud region IILoren et aLllTMOh . 



ray, UV, and accretion heating. iRoberge et al.l ()2001[ ) 
detected H2 in absorption towards AB Aur, probably 
located in the envelope around the star. They derive a 
temperature and column density of T = 212 K and N(H2) 
= 6.8 X 10^^ cm^^ for the absorbing gas. If we assume 
the emission is spread evenly over our TEXES beam, the 
flux in our lines due to this gas is less than 10~^^ ergs 
s~^ cm~^, sig nificantly lower than our detected fluxes. 
iBrittain et al.l ([2003) detected CO fundamental rovibra- 
tional emission from AB Aur and found that the emission 
was coming from both hot (1540 K) CO in the inner disk 
rim and cold (70 K) gas farther out in the flared region 
of the disk. Our derived gas temperature based on the 
II2 lines falls between the hot and cold components seen 
in the CO observations, suggesting we are not seeing the 
same gas. 

4.1.2. DoAr 21 

DoAr 21 is located in the p Ophiuchus molecular 
cloud and is classified as a weak- line T Tauri star 
(|Bouvier fc Appenzellei|[l993 ). The distance to p Ophi- 
uchus remains a source of debate wit h distance csti- 
mates ranging b etween 119 and 165 pc (fMamaieljl2008l : 
iLombardi erall 120081 : iLoinard et all l2008f) . Fbr the 
sources in our sample, we a dopt t he recent distance es- 
timate by ILombardi et al.l ()2008l ) of 119 pc. Contin- 
uum observations r eveal only a slight infrared excess 
(jWilking et al.|[200ll ) and millimeter observations suggest 
the amount of circumstellar dust is less than 10~^ Mq 
corresponding to a gas mass less t han 10~^ Mfn assum- 
ing t he standard gas-to-dust ratio ([Andre fc Montmerld 
1994f). Anv circum stellar disk around the star is tenuous. 



Barv et al l (|2002l ) detected narrow, 9 km s'^ FWHM, 



near-infrared H2 v=l-Q S(l) emission centered at the 
systemic velocity of DoAr21 and concluded that the line 
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arose from a circunistellar disk within ^^110 AU of the 
star. Due to the a bsence of a double-peaked hue profile, 
iBarv et aTl (|2002D concluded that the circumstellar disk 
in the system has an inclination <45°. iBarv et all (120031 ) 
calculated that the line emission could be produced by 
a mass of 2.7 x 10"^ Ha gas in LTE at T = 1500 
K located in a thin surface layer of the disk. This mass 
was calculated assuming a distance of 160 pc. Correcting 
their mass to the 119 pc we assume for DoAr 21, this gas 
would produce a S(2) hne flux at 12 ^m of 0.6 x 10~^^ 
ergs s~^ cm^^, smaller than both our detection from the 
IRTF and the upper limit based on Gemini data. 

We observed DoAr 21 with TEXES on both the IRTF 
and Gemini North telescopes. From the IRTF, we de- 
tected narrow, FWHM ^^6 km s^^, II2 S(2) emission 
(Figure E]), while the line was not seen in our Gemini 
observations. Since our slit size on the sky is smaller on 
Gemini than the IRTF, a plausible explanation for the 
discrepancy is that the emission is spatially outside our 
Gemini slit. At the II2 S(2) setting, our slit widths are 
1.4" and 0.54" on the IRTF and Gemini respectively. 
At an adopted distance of 119 pc to DoAr 21, these 
slit widths translate to physical distances from the cen- 
tral source of 83 AU and 32 AU. Our data suggest that 
the observed H2 emission arises outside of the inner 32 
AU around DoAr 21. The ro vibrational emission seen 
in DoAr 21 (jBarv et al.l 1200^ was obse rved using the 
Phoenix spectrometer llHin kle et al.ll 19981 ) on the NOAO 
4-meter telescope on Kitt Peak with a seeing-limited slit 
width, the same as TEXES at the IRTF. If our obser- 
vations probe the same gas, this is consistent with the 
possibility that the emitting region is outside of the in- 
ner ~30 AU. 

4.1.3. Elias 29 

Elias 29 is a class I protostar located in the p Ophi- 
uchus molecular cloud. Along with several nearby proto- 
stars, it is located along a dense ridgelike stru cture seen 
in HCO+ 3-2 emission ()Boogert et al.ll2002aD . The cen- 
tral source is surrounded by a circumstellar disk and a 
remn ant envelope several times more massive than the 
disk (jLommen et al.ll2008l). Modehng of the spectral en- 
ergy distribution (iBoogert et ahl l2002aD constrains the 
disk size to ^^500 AU with an inclination less than 60° 
and mass of 0.012 Mq. 

Using the Infrared Spac e Observatory ( I SO) Short 
Wavelength Spectrometer, ICeccarelli et al.l (|2002[ ) saw 
evidence for a disk around Elias 29 with a superheated 
surface layer and mass similar to disks around Her- 
big AeBe stars. These authors suggested that Elias 29 
may actually b e a de eply embedded Herbig AeBe star. 
IBoogert et al.l (|2002al ) also suggested Elias 29 might be 
a heavily extincted T Tauri or H erbig AeBe star. Elias 
29 drives a bipolar CO outflow (|Bontemps et al.l 119961) 
with v elocities approaching ^ 80 km s~^ (jBoogert et al.l 
l200 2aD. Knotty H2 1-0 S(l) emission suggests the 
presence of a processing jet interacting with the sur- 
rounding medium a nd clearing the protostellar envelope 
(jYbarra et al.l[2006l ). 

We used TEXES on the IRTF in 2003 June and on 
Gemini North in 2006 July to observe Elias 29 at the 
H2 S(l) and S(2) settings. Wc obtained observations of 
the H2 S(9) hne taken with Kcck/NIRSPEC during sev- 
eral runs between 2000 and 2005. The data are shown 
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Fig. 6.— Elias 29 NIRSPEC/Kcck S(9), TEXES/IRTF S(2), 
S(l), and TEXES/Gemini S(2), S(l) overplotted with 1-component 
mo del fit. The dashed lin e indicates the systemic velocity of Elias 
29 l IBooeert et al.l l2002al ). The continuum has been scaled to 
agree with the ISO SWS measurement. The bump redward of 
the S{9) feature is due to telluric contamination. Poor telluric di- 
vision of a water feature is apparent on the blue side of the Gemini 
S{2) line. 

overplotted with a single temperature LTE model fit in 
Figure [SI The Gemini spectrum shows poor telluric di- 
vision on the blue side of the S(2) line. A telluric water 
feature at this velocity apparently was not well corrected. 
The results of the single temperature LTE model fits to 
the data are listed in Table [31 The data are well fit by 
emission from less than 1 of gas at T'^IOOO K. The 
S(2) line was >3 times stronger in our IRTF observa- 
tions compared to Gemini and the S(l) line is clearly 
detected in our IRTF observations but not seen from 
Gemini. There are two possible explanations for these 
discrepancies. Either the H2 line flux actually changed 
between our observations from the IRTF and Gemini, or 
the emission is spatially extended and our wider slit on 
the IRTF took in more of the line flux. We saw no evi- 
dence of spatially resolved H2 emission along the slit of 
our Gemini observations at scales of 0.4 arcsec (^^25 AU 
in radius at 119 pc). The continuum level did not vary 
significantly between our IRTF and Gemini observations. 

4.1.4. GSS 30 IRS 1 

CSS 30 IRS 1 is a class I source (IWilking et al.l[l98l 
in the p Oph iuchus molecular cloud. Ne ar-infrared po- 
larimetry by IChrvsostomou et"aLl (|1997') suggests that 
the source is surrounded by a large flattened dusty en- 
velope and a more compact circumstellar disk. Model- 
ing of near-infrared polarimetry data suggests the disk 
has an inclination of ^6 5°, i.e. closer to edge-on 
(jChrvsostomou et al.lll996[ ). Nearby molecular outflow 
activity has been observe d but is not c l early associated 
with GSS 30 IRSl itself. iTamura et all (|l990h observed 
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Fig. 7.— GSS 30 NIRSPEC/Keck S(9) and TEXES/Gemini 
S(4), S(2), S{1) overplotted with 2-component model fit. The 
dashed line shows the systemic velocity £]Pontoppidan et al. 2003)- 
The continuum has been scaled to agree with the Spitzer IRS mea- 
surement. 

high velocity miUimeter CO emission to the south of GSS 
30 IRSl that is hkely asso ciated with the nearby VLA 
1 623 let (lAndre et al.lll990l) 

iPon toDpidan et aLl ( 2003 ) described observations of 
fundamental rovibrational CO emission at 4.7 /im. The 
lines are unresolved at R = 5000 and are spatially ex- 
tended up to 320 AU from the central source. The au- 
thors proposed that the line emission arises from post- 
shocked gas from the inner region of the circumstellar 
disk, which is then scattered into our line of sight by the 
surrounding envelope. The lines are consistent with 1- 
100 Mg of gas in LTE at T = 515 K with a spatial extent 
of -20-100 AU. 

Our pure rotational H2 data cannot be fit by a LTE 
single temperature/mass model, but is fit reasonably 
well with a two component model (Figure [7]). Using a 
two component fit, we f ind a temperature for the low J 
lines nearly identical to iPontoppidan et al.1 (|2002D . The 
S(9) emission arises in a small amount of hot gas, sig- 
nificantly hotter at 3300 K than the dust evaporation 
temperature, and has a broader line width than the 
other H2 lines. Because such hot gas temperatures are 
only effectively constrained by the S{9) line, our deter- 
mination of this temperature is more uncertain. Since 
the H2 molecule may form in a high rotational level 
([Wagenblastl 1199^ . this apparent high temperature, if 
va lid, may be due to II2 for mation pumping. As noted 
bv lPontoppidan et al.l ()2002f ) , the intrinsic line widths are 
expected to be ~10 km s~^ or larger in the case where a 
dissociating accretion shock is responsible for the emis- 
sion. The width of the low J pure rotational H2 lines 
are consistent with this model. The hot gas, however, is 
likely located in a separate location. 




50 
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Fig. 8.— GV Tau N TEXES/Gemini S(2) overplotted with 
Gaussian fit. The d ashed line indicates the systemic velocity of 
the GV Tau system IIHogerheiide et aLllTOOa) . 



4.1.5. GV Tau N 

GV Tau is a pre-main-sequence binary system in the 
L1524 molecular cloud. The two com ponents are sep- 
arate d by 1.2" (170 AU at 140 pc) (|Leinert fc Haas! 
Il989t ). The southern component is optically visible while 
GV Tau N is heavily extincted. Near-infrared imaging 
and optical polarimetry show that the system is sur- 
rounded by a flattened, edge-on circumbina ry envelope 
out to -1000-1500 AU (iMenard et al.lfl99l . GV Tau 
N shows near-infrared variability on timescales as short 
as a month, which has been attributed to dum piness 
in the surrounding material ( [Leinert et al.|[200ll ). The 
GV Tau system has a spectral energy distribution ris- 
ing through the mid-infrared leading to its classification 
as a class I source , howe ver, millimeter observations by 
iHogerheiide et al.] (|1998D suggest that most of the enve- 
lope has disappeared. GV Tau N appears to be driving 
a Her big-Haro outflow (jPevine et al.lll999f) . iGibb et all 
()2007l ) reported the detection of near-infrared absorption 
hues due to CO, HCN, and C2H2 toward GV Tau S and 
derived a CO rotational temperature of —200 K, suggest- 
ing the o bserved gas is in t he inner region of a circumstel- 
lar disk. iDoppmann et all ([2008) observed HCN, C2II2, 
and CO absorption towards GV Tau N but did not detect 
molecular absorption towards GV Tau S. Mid- infrared 
observations with TEXES on Gemini in 2006 November 
also showed HCN absorption towards GV Tau N but not 
GV Tau S (Najita et al., in preparation). 

We observed GV Tau N at the S(2) setting (Figure [5]) 
from Gemini and detected emission at a flux level and 
line width consistent with other detections in our sam- 
ple. The detection of just a single line precludes an 
esti mate of the tempe rature and mass of the emitting 
gas. IDoppmann et al] ([2008i) observed near-infrared H2 
emission from GV Tau N in the v—1-0 S(2) and v—1- 
S(0) lines. The line centroids are consistent with our 
mid-infared S(2) detection and with the systemic veloc- 
ity of GV Tau. The u=l-0 S(2) line emission is stronger 
than the v—1-0 S(0) by a factor that is consistent with 
the shock excitation seen in classical T Tauri stars by 
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Fig. 9.— HL Tau NIRSPEC/Keck S(9) and TEXES/Gemini 
S(4), S(2), S(l) overplotted with 2-componcnt fit. The dashed line 
shows the stellar velocity iWhite £: H illenbrand 2004). The con- 
tinuum has been scaled to agree with the ISO SWS measurement. 
The bump redward of the very weak S(9) feature is due to telluric 
contamination. 



iBeck et al.l (l2008l) . 



4.1.6. HL Tau 



The young stellar object, HL Tau, is in an interme- 
diate stage between an embedde d class I protost ar and 
an optically visible T Tauri star (|Pvo et al.ll2006f ). It is 
surrounded by a 0.05-0.07 Mq circumstellar disk with 
an outer radius of 90-160 AU (iMundv et al.lll996f) in - 
clined by 66-71 degrees from face-on ( Lucas et al.ll2Q04l ). 
Molecular carbon absorption suggests the presence of an 
infalling envelope (jGrasdalen et al.l ll989'). A coUimated 
jet is seen from HL Tau in [Fe H] 1.6 /xm emission sur- 
rounded by a wide-angled wind, which produc es shocked 
near-infrared Hq emission (Takami et al 2007^. Spatially 
resolved observations by Beck et al. ( 2008) with NIFS on 
Gemini North confirm that the near-infrared H2 v=l-0 
S(l) emission is consistent with the location of the jet 
associated with the system. Near ly half of the NIR H2 
emission seen bv'Beck e t al.l ()2008 l) is spatially coincident 
with the continuum. If the spatial distribution of the 
mid-infrared H2 emission is the same as the near-infrared 
H2 emission, our observed emission lines may also origi- 
nate in circumstellar gas shocked by the jet. Broad near- 
infrared CO emission from the hot (T ^ 1500 K) inner 
disk along with narrow CO absorption likely originat- 
ing in the outer flar ed disk are obser v ed fro m HL Tau 
(jBrittain et al.ll2005f ). iBrittain et all (|2005f l estimated 
that the emission arises from CO with a column density 
of 4 X 10^^ cm~^ in a region between 0.066 AU and 0.53 
AU. This translates to ~ lO"* of H2 at 1500 K, which 
would produce pure rotational mid-infrared H2 emission 
at levels ~ 10^^^ ergs s^^ cm^^, well below our detec- 
tion limits. This combined with the average FWHM of 



the CO lines observed by IBrittain et"al1 ()2005[ ) (45 km 
s~^) makes clear that our observations are probing gas 
at larger radii. From the IRTF, we detected emission 
in the H2 S(2) line, while from Gemini we detected H2 
S(l), S(2), and S(4) emission. A single temperature LTE 
model comprised of ~1 of 465 K gas fits the observa- 
tions well (Figure O. The emission lines are all narrow 
with FWHM near 10 km s"^ 

4.2. Putting TEXES Results into the Context of Other 
Circumstellar Cas Observations 

Emission has been detected in the near- infrared ?;=1- 
S(l) transition of H2 from several T Tauri stars 
(|Barv et al.ll2003l ). The detected emission shares similar 
characteristics with our observations of the mid-infrared 
H2 lines. In both cases the lines are narrow and centered 
at the systemic velocity of the star. The observed near- 
infrared H2 lines all have FWHM ~10 km , suggesting 
that the emission arises 10-50 AU away from the star and 
is possibly the result of excitation by UV or X-rays. The 
spatially resolved observations of v—l-Q S(l) H2 emis- 
sion from the circumstellar env ironments of classical T 
Tauri stars bv lBeck et al.l ()2008l ) are most consistent with 
shock excited emission from outflows rather than UV or 
X-ray excitation. 

Three of the sources with dete cted near-infrared H2 
emission in the iBarv et al] (|2003f l sample were also ob- 
served as part of our TEXES H2 mid-infrared survey. 
Two of those sources, GG Tau and LkCa 15, do not 
have detectable levels of mid-infrared H2 emission while 
the other, DoAr 21, was detected in the S(2) line in 
our IRT F observations bu t not when re-observed from 
Gemini. iBarv et al] ()2003D computed the amount of H2 
gas required to produce the observed emission in these 
sources under the assumption of LTE at T=1500 K. 
Based on their derived gas masses for these three sources, 
we computed the amount of emission we would see in the 
mid-infrared H2 S(l) and S(2) hues. For GG Tau and 
LkCa 15, the predicted line fluxes are all less than a few 
times 10~^^ ergs s~^ cm~^ consistent with our derived 
upper limits. In the case of DoAr 21, the predicted line 
fluxes are less than 10~^^ ergs s~^ cm~^, smaller than 
both our non-detections from Gemini and the detected 
S(2) line flux observed with the IRTF from DoAr 21. 

CO fundamental rovibrational emission near 4.6 //m 
is detected in many T Tauri stars (Naiita et al. 200^ 
and Herbig AeBe stars (Bl ake fc Boogertll2004D . The A- 
values for the CO fundamental lines are much larger than 
those of the pure rotational mid-infrared H2 lines and so, 
assuming LTE, the CO lines are more sensitive to small 
column densities of gas. The lines are broad and centrally 
peaked with FWHM of 50-100 km s~^, s uggesting the 
emissi on arises from < 0.1 AU to 1-2 AU (jNaiita et al.l 
l2007al) . The CO near-infrared observations typically re- 
veal temperatures of 1000-1500 K and C O column densi- 
ties of ~10^^ cm~2 (jNaiita et al.ll2007al) . Assuming the 
CO fundamental emission arises from 0.1-2 AU with a 
CO column de nsity of 10^^ cm ~^ and a CO/H2 ratio 
of 2.7 X 10""* (iLacv et al.lll994l ) gives a mass of H2 in 
this region of ~10~^ M^. The line fluxes in the mid- 
infrared H2 lines from such gas at T = 1000 K are less 
than 10^^^ ergs s~^ cm~^, smaller than our detection 
limits. The additional line broadening due to the larger 
rotation speeds in this part of the disk would further 
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decrease the chances of this gas being seen in our high 
resolution mid-infrared H2 observations. HD 141569 is a 
unique source among those with observed near-infrared 
CO emission in that the emission arises from gas at a 
much cooler temperature (1 90 K) and a t a location > 
17 AU in the disk (Brittai n et alil2003f ). iBrittain et al] 
(|2003f) derived a mass for the emitting CO gas of 10 
g (~10~^ M0). For a rotational tem perature of 190 K , 
assuming the CO/H2 ratio derived bv lLacv et all (|1994j ). 
the expected mid-infrared II2 line fluxes are smaller than 
10~^° ergs s~^ cm~^, consistent with our non-detections. 

4.3. Location of Emitting H2 Gas and Possible 
Excitation Mechanism 

The six stars with detected H2 emission all have nar- 
row, spectrally resolved line widths between 7 and 15 km 
s~^. Added to the fact that the line fluxes are all similar 
and that the lines are centered at the stellar velocity, this 
suggests that the excitation mechanism is similar in each 
case. If the emission originates in a circumstellar disk, 
our spectrally resolved lines allow for the determination 
of the approximate emission radius. We created simple 
models for the line widths originating from a Keplerian 
disk that contributes equally at all radii within some an- 
nulus and convolved these line profiles with the TEXES 
instrumental profile and thermal broadening appropriate 
for gas at T = 500 K, roughly the temperature derived 
from our observations. For a star of one solar mass with 
a disk inclination angle of 45°, the range of line widths 
detected in our sample corresponds to emission from disk 
radii between 10 and 50 AU. 

In all cases, our observations show the emission is spa- 
tially unresolved along the slit of our Gemini observations 
at scales of 0.4 arcsec and coincident with the source 
continuum. This combined with the fact that the lines 
are spectrally resolved with FWHM ~10 km s~^ sug- 
gests that we are not seeing emission from an extended 
envelope surrounding the source, which would produce 
narrow, spectrally unresolved lines. Furthermore, it is 
unlikely we are seeing the results of shocks associated 
with jets and outflows commonly seen near young stars, 
as that would produce broader lines with some displace- 
ment from the stellar velocity. 

However, the spatially resolved (~0.1") observations of 
near-infrared H2 emission from the circumstellar environ- 
ments of six classical T Tauri stars by Beck et al. (2008) 
demonstrate that we cannot rule out shocked emission 
as the source of our observed mid-infrared lines in all 
cases. One source we sha re with lBeck et al.l ()2008[) is HL 
Tau. iBeck et all (|2008f) concluded that their spatially 
resolved observations of the near-infrared v—1-0 S(l) H2 
line from HL Tau are consistent with the location of the 
jet in the system. In addition, nearly half of the ob- 
served H2 emission from HL Tau is spatially coincident 
with the continuum and the line centroid is within 10 
km s~^ of the stellar velocity. Assuming the spatial po- 
sition of the near-infrared and mid-infrared H2 emission 
is the same, this suggests that the emission we observed 
may also arise in shock excited circumstellar gas. How- 
ever, there is a notable difference in the overall results 
of the spatiall y resolved, near-infrared H2 observations 
by IBeck et al.l (2008) and our mid-infrared H2 observa- 
tions. In half of the six stars observed, the near-infrared 
line centroids differed by more than 10 km s^^ from the 



stellar velocity whereas in our sample, all of the H2 line 
centroids are within a few km s^^ of the stellar velocity. 

The results of spectral energy distribution modeling 
of the dust temperature in the surface layer of a disk 
around a typical T Tauri star show that the dust tem- 
perature is < 200 K at dis k radii larger than 10 AU 
where our H2 emission arises (iChiang fc Goldreichll997b . 
Stellar heating of dust grains in the disk atmosphere 
coupled to the gas temperature through gas/grain colli- 
sions is insufficient to explain th e high gas temp e ratures 
derived from our observations. iCarmona et al.l (|2008ar ) 
computed the expected emission in the H2 S(l) and 
S( 2) lines from the opt i cally t hick, two-layer disk model 
of IChiang fc GoldreichI (|1997T ) in which dust in the disk 
surface layer absorbs stellar radiation and heats the gas. 
They found that the amount of gas in the warm sur- 
face layer of the disk is less than a few Earth masses. 
At a distance of 140 pc, this leads to predicted mid-IR 
H2 line fluxes of 10^^^ — 10~^^ ergs s~^ cm^^, much 
lower than bot h the leve ls of our detections and our up- 
per limits. C armona et a l. (2008a) pointed out that the 
two-layer approximation to the disk structure leaves out 
details that could significantly contribute to H2 emission. 
They found that departures from thermal coupling be- 
tween gas and dust in the disk surface layer as well as 
larger than interstellar gas-to-dust ratios can lead to de- 
tectable levels of H2 emission. Two plausible mechanisms 
for additional gas heating in the surface layers of disks 
are accretion shocks due to infalling matter onto the disk 
and X-ray/UV irradiation. That accretion onto the disk 
may play a role in exciting H2 emission is consistent with 
the preferential detection of H2 emission from the class 
I sources in our sample, which possess a surrounding en- 
v elope of material in addition to a disk. 

iNeufeld fc HollenbachI (|1994f ) have calculated the phys- 
ical and chemical structure of shocks resulting from ac- 
cretion onto a circumstellar disk, but they did not pub- 
lish the strength of the pure rotational H2 lines from 
their model. However, we have reviewed the results of 
unpublished models from their study to determine the 
parameters that would produce H2 emission at the lev- 
els de tected in our observations. iNeufeld fc HollenbachI 
()1994[ ) only considered preshock densities above 3 x 10^ 
cm~^. However, below this density, the fraction of the 
total cooling from H2 emission increases. A shock with 
a preshock density of 10^ cm~^ striking the disk at 5 
km s^^ at 30 AU implies an accretion rate of 5 x 10~^ 
Mq yr"^ The resuhing emission in the S(l), S(2), and 
S(4) H2 lines is ^^10^^ Lq. This accretion rate is similar 
to the measured rates for the sources where we have de- 
tected H2 emission and these line luminosities are only 
slightly lower than our observations. One way to increase 
the H2 cooling in the model is to suppress the coohng in 
H2O lines in the shock by assuming that H2O is frozen 
out as water ice in the preshock gas. It is plausible that 
the H2O would be frozen out for the 5 km s~^ shocks 
since they are slow enough that they would not return 
the H2O to the gas phase. If this were the case, the 
luminosity in our H2 lines would increase by an order 
of magnitude to about 10~^ Lq. The line luminosities 
detected in our sample range from 10~^ to 10~^ L0 sug- 
gesting that shock heating due to accretion onto a disk 
is a plausible excitation mechanism. Additional details 
including predicted CO line strengths can be found in 
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Fig. 10. — Gemini H2 line luminosities vs. X-ray luminosity of 
each source. The solid points are H2 detections and the arrows 
represent upper limits. No correlation between the presence of H2 
emission and X-ray luminosity is apparent. 

t he Appendix. 

iNomura et al.l (120071 ) modeled the molecular hydrogen 
emission from a disk surrounding a typical T Tauri star, 
taking into account the heating of gas by X-ray and UV 
irradiation from the central star. The resulting gas tem- 
perature in the surface layer of the disk is much higher 
than the dust temperature. X-ray heating dominates in 
the inner region and surface layer of the disk. At 10 
AU, the disk surface temperature reaches over 1000 K. 
Even at 100 AU, the temperature in the disk surface layer 
reaches 200 K. These temperatures combined with the 
fact that the gas is hotter than the dust creates favorable 
conditions for the production of mid-infrar ed H2 emission 
lines. The predicted line fluxes derived bv lNomura et al.l 
(|2007f ) appropriate for TW Hya vary depending on the 
adopted dust size distribution. The quoted line fluxes 
at our wavelengths are 4.4 — 9.9 x 10~^^ ergs s~^ cm~^ 
at S(l), 1.7 - 4.8 X 10-15 g^-gg g-i ^^-2 ^(^2), and 
0.4 - 8.8 X 10-15 ergs s-^ cm-^ at S(4). At distances 
more typical of the sources in our sample, these fluxes are 
~10-i^ ergs s-i cm-^. In the sources where we detect 
line emission, the derived line fluxe s are all larger, typ 



ically 10-15-10-" ergs s'^ cm-^. iGorti fc HollenbachI 
(|2008l ) have also modeled line emission from the upper 
layers of optically thick disks. Their predicted H2 emis- 
sion line luminosities are 10-^ to 10-^ Lq and are consis- 
tent with the range of detected emission in our sample. 
The difference in the predictions of the two models shows 
that the modeling is very sensitive to the input param- 
eters. The FUV flux assumed by iGorti fc HollenbachI 
()2008f l in their standard model is five times higher and 
th e X-ray flux is three t imes larger than the values used 
by iNomura et al.l ()2007f ) which leads to the higher pre- 
dicted H2 line fluxes. 

In Figure [TOl we present a plot of our H2 detections 
and upper limits versus the X-ray luminosity of the stars 
in our sample. If high X-ray luminosities were responsi- 
ble for heating the gas in disks above the dust temper- 
ature to produce detectable II2 emission, we should see 
a correlation between detected II2 emission and X-ray 



log FUV Luminosity (L^) 

Fig. 11. — Gemini H2 line luminosities vs. FUV luminosity of 
each source. The solid points are H2 detections and the arrows 
represent upper limits. There is a hint of a correlation between 
FUV luminosity and detected H2 emission for the small sample of 
S{4) data points, however there is no clear overall correlation. 
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Fig. 12. — Gemini H2 line luminosities vs. accretion rate mea- 
sured in each source. The solid points are H2 detections and the 
arrows represent upper limits. The sources with detected H2 emis- 
sion fall in the middle of the range of accretion rates. 

luminosity. As seen in Figure llOi such a correlation is 
not apparent in our data. In several cases, we detect II2 
emission from sources at the faint end of the distribution 
of X-ray luminosities, while we did not detect II2 emis- 
sion from the most X-ray luminous stars in our sample. 
It is worth noting that X-ray flux has been observed to 
vary by more than an order of magnitude on timescales 
of days ( Kastner et al.|[l"999f ) . To definitively test for ob- 
servational evidence of a correlation between X-ray flux 
and II2 emission requires a series of coordinated obser- 
vations. UV irradiation of the disk from the central star 
may also contribute to heating the gas in the disk, leading 
to detectable levels of II2 emission. We plot II2 line lu- 
minosity versus FUV luminosity in Figure [TT] and find no 
clear correlation, however. Since accretion heating may 
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play a role in producing H2 emission, we plot H2 line 
luminosity versus accretion rate in Figure [12] to look for 
a correlation. The sources in our sample with detected 
H2 emission are located in the middle of the range of 
accretion rates, showing no clear correlation. Since the 
dynamic range separating our H2 detections from the 
non-detections is not very large, the lack of a correlation 
between Lx, Lpjjv, M and detected H2 does not con- 
clusively rule out the possible importance of X-ray, UV, 
and accretion heating in producing H2 emission. We also 
searched for a correlation between H2 emission and mass, 
age, and inclination but found none. 

5. CONCLUSIONS 

We have carried out a survey for pure rotational H2 
emission from the circumstellar environments surround- 
ing a sample of 29 stars with disks and detected emission 
from 6. In the case of non-detections, our upper limits 
constrain the amount of T > 500 K gas in the surface 
layers of the circumstellar disks to be less than a few 
Earth masses. Several objects in our survey have tran- 
sition object SEDs implying the presence of an optically 
thin, dust-depleted inner disk: DoAr 21, GM Aur, HD 
141569, and LkCa 15. Among these sources, only DoAr 
21 shows H2 emission and it appears to be far from the 
star. One pos sible explanation f or transitional SEDs is 
gram growth (jStrom et al.lll989t iDullemond fc DominikI 
|2005() whereby the inner disk becomes optically thin yet 
remains gas rich. This gas could be heated through accre- 
tion as well as X-ray and UV heating. GM Aur is of par- 
ticular interest since it has an accretion rate (^10^* M© 
yr~^) similar to typical CTTS accretion rates so should 
have as much gas in the inner disk as a typical CTTS. 
If grain growth is responsible for the transitional SEDs 
of these sources, the available heating mechanisms are 
insufficient to produce detectable H2 line emission. Al- 
ternatively, grain growth may not be a good explanation 
for a transitio n object SEP, as sug gested by other demo- 
graphic data (iNaiita et al.ll2007bf ). 

In all cases, the detected emission lines are narrow and 
centered at the stellar velocity. The narrow range of line 
widths, FWHM between 7 and 15 km s~^, along with the 
fact that the line fluxes are all similar, suggests that the 
mechanism for exciting the emission may be the same in 
each case. Four of the six targets with detected emission 
are class I sources that show evidence for surrounding 
material in an envelope in addition to a circumstellar 
disk. It is possible, and likely in the case of HL Tau, 
that the H2 emission we observe is a result of gas in the 
circumstellar envelope being shock heated by an outflow. 
However, the fact that all of the H2 line centroids in our 
sample are within a few km of the stellar velocity 
argues against this being the case for all of our detections. 

Under the assumption of emission from a disk in Keple- 
rian rotation, the narrow line widths imply that the emis- 
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sion arises at disk radii from 10-50 AU. At such large disk 
radii, additional heating of the gas besides heating due 
to collisions with dust grains is required to explain the 
temperatures derived from our II2 observations. Both 
X-ray/UV irradiation of the disk surface layer and accre- 
tion shocks resulting from matter infall onto the disk are 
plausible candidates. With the exception of DoAr 21, 
all of the sources where we detect H2 emission possess 
both a circumstellar disk and a surrounding envelope of 
material. This lends support to the possibility that the 
II2 emission we observed may be the result of shocks in 
the disk due to infalling material. 

Models of molecular hydrogen emission from disks that 
assu me sufficient levels of stella r X-ray and UV irradia- 
tion (|Gorti &: HollenbachI 120081) predict line fluxes that 
are consistent with our observations. In contrast, mod- 
els which assu me smaller values of stellar UV and X- 
ray irradiation (|Nomura et al.|[2b07f) produce weaker II2 
emission than observed in our sample. We looked for 
evidence of a correlation between X-ray/UV luminosity 
and the presence of H2 emission but found none. We note 
that the X-ray and UV luminosities used for the purpose 
of searching for a correlation with II2 emission were not 
measured at the same time. To definitively test for a 
correlation between X-ray/UV luminosity and the pres- 
ence of II2 emission will require a series of coordinated 
observations. 
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Physical properties of the stars in our sample 
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15''49™57^ 


.7 


-03°55'17".0 


B9.5 


HAeBe 


994 


> 10.0^ 


< 28.132 


0.4863 


0.0438 


HD 163296 


17'^56™21^ 


.3 


-21°57'23".0 


Al 


HAeBe 


1224 


6.02 


29.632 


0.2673 


0.00133 


HL Tau 


04'"31™38^ 


.5 


18°13'58".0 


K9 


I 


14024 


0.7734 


30.61° 


0.17 


1.629 


IRAS 04278+2253 


04^30™ 50^ 


.7 


23°00'll".l 


Fl 


I 


14024 






7.1 


66.129 


L1551 IRS 5 


04'*31™34^ 


.2 


18°08'05".3 


G-K^s 


I 


14024 




28.310 




< 14036 


LkCa 15 


04*^39™ 17^ 


.8 


22°21'03".5 


K5 


III 


14024 


11.72 


< 22.625 


0.005637.38 


0.06839 


Lk Ha 225 


20'"20™30^ 


.8 


41°21'24".8 




HAeBe 


10004" 










MWC 758 


05'"30™27^ 


.4 


25°19'56".8 


A3 


HAeBe 


200^ 






0.0733 




RW Aur 


05/107^493 


.6 


30°24'05".4 


G5 


I 


14024 


2.5734 


< 29.441 


0.0653.13 


15.816 


SVS 13 


03'"29™03^ 


.6 


31°16'01".2 




I 


30042 




< 30.343 


0.67 


6.344 


TW Hya 


ll''01™51^ 


.9 


-34°42'18".3 


K7 


II 


5145 


1046 


29.747 


0.0017 


0.01848.49 


V892 Tau 


04'*18™40^ 


.7 


28°19'16".2 


A6 


HAeBe 


14024 




3O.9IO 


0.36 




V1057 Cyg 


20'"58™53^ 


.1 


44°15'28".6 




fuori 


60021 




< 31.017 


0.213.23 




V1331 Cyg 


2l'"01™09^ 


.1 


50°21'45".2 


G5 


fuori 


55050 


0.823 


< 3O.9I2 






VV Ser 


18''28™47^ 


.8 


00°08'39".7 


A2 


HAeBe 


24551 






5.123.8 


4.578 


Z CMa 


07/io3'"42« 


.0 


-11°33'02".8 


F69 


fuori 


115052 




30.332 


5.963.53 


79021 



References. — l |Javawardhana et all 1120011' ). TThi et al.l 1 1200 ID . Walenti et all ll2003l ').4Hipparcos. 1van den Ancker et al.l 

||200T),6Hcrbcrtz ct al.' fT99T),7BcrEhofcr el's] ||T996| ').^ Garcia Lopez ct al.' (2006), ^van den Ancker ct al. ( 1998),iOGiidcl cTaO 
^007),iiLombardi ct aL (2008 ,), 12 'Walter fcTCuhil (19 81),13 Valcnti ct al. ( 1 993), i4Tokunaga ct al, ( 2004 ),15 Prate ct aQ 
(2003),i''HartiKan ct al] llT995l) .i7|Carkncr c^g__g99i) , i* Luhman & Rickc l ll9^ .i'^Imanish i et al.l (2002),20Natta ct al] 
f2006 ).2iHartniann fc Kenyon! |T996), 2= Skinner et al.l l'^006),23Hcrbig & DaliinI pb06),24Elia3 jl978),25Ncuhauscr ct al] 
llT995l ),2e|G-ullbring ct al. ( 1998),27FciKclson & DcCampli ( 1981),2*^ Gagnc ct al. (2004),2SWhitc & Hillenbrand (2004),3':Dolan & Mathici] 
Il2p0l!),3l|c5;ivct ct al. (2004),32Stclzcr ct al. (2006b), 33Swartz ct al. ( 2 00^,34Siess ct a l. ( 1999).3 5,Mundt et al.. tl98 5').36Fullcr ct af] 
llT995),37Bcrgin ct al. (2004),38Espaillat ct at. ( 20071.3^ Hartmann eTaTI ||T998| ), 4'- Marvel ll2005l) ,4 l|Gahml ( 1980),42 Ccrnii 
119901 ) ,4^tman et al. (,2002l),4lEdwards et al.l l l2003i),4 5 iMamaieia ||2005|' ). '^\\}c\a.A& et al.l ll200l),4'TStelzer fc Schmit^ 
ll200^ ),4Wencar & Batalhi ll25ol ).4' Blerczeg et al.l ll200i ).5' 15hevchenko et"all llT99ll ').5l [Uhavarria et all 119881' ). 5 TOerbst et al.l 
lT978l) .53 ;gtelzer et al.i (l2006a^ 

If not otherwise noted, value derived using TW Hya lUE Spectrum then scaling by relative accretion rates. 
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TABLE 2 
Summary of Line Detections 



Line Equivalent 



Star 


Instrument 


A 
(Mm) 


Date 


Gont. 

(Jy) 


Line Flux'' 


Luminosity 
(10-6 Lq) 


Width 
(km s-i) 


FWHM 
(km s-i) 


AB Aur 


TEXES /IRTF 


12.279 


dec02,dec03 


12.7'' 


0.93 (0.25) 


6.0 


0.86 


(0.23) 




7.0 


AB Aur 


TEXES/Gemini 


8.025 


nov06 


12. 7*" 


1.47 (0.34) 


9.5 


0.93 


(0.21) 




10.4 






12.279 


nov06 


14. 7*" 


0.53 (0.07) 


3.4 


0.44 


(0.06) 




8.5 






17.035 


nov06 


24.6*' 


1.10 (0.3) 


7.1 


0.76 


(0.21) 




9.0 


AB Aur 


TEXES/Gemini 


8.025 


oct07 


12.7'' 


< 1.24= 


< 8.0 














12.279 


oct07 


14.7^ 


0.56 (0.07) 


3.6 


0.47 


(0.06) 




9.4 






17.035 


oct07 


24.6'' 


0.57 (0.16) 


3.7 


0.29 


(0.08) 




6.5 


AB Aur 


NIRSPEC/Kcck 


4.695 


jan01-dcc02 


9.6'' 


0.90 (0.07) 


5.8 


1 02 








DoAr 21 


TEXES/IRTF 


12.279 


jun03 


0.14'' 


0.33 (0.09) 


1.5 


24.06 


(6.46) 




5.6 


Elias 29 


TEXES/IRTF 


12.279 


jun03 


24.0^= 


2.45 (0.33) 


10.9 


1 97 


(C\ 1 7^ 
[\J. 11) 




10. 1 






17.035 


jun03 


30.5^= 


1.64 (0.30) 


7.3 


0.94 


(0.18) 




11.9 


Elias 29 


TEXES/Gemini 


12.279 


juloa 


24.0^= 


0.70 (0.12) 


3.1 


0.36 


(0.06) 




12.7 


Elias 29 


NIRSPEC/Kcck 


4.695 


julOO-aprOS 


22.9'' 


2.12 (0.33) 


9.4 


0.44 


(0.07) 




21.8 


GSS 30 


TEXES/IRTF 


12.279 


jun03 


14.7^* 


1.19 (0.25) 


5.3 


0.98 


(0.21) 




6.8 


GSS 30 


TEXES/Gemini 


8.025 


juloa 


4.21^* 


0.78 (0.07) 


3.5 


1.50 


(0.14) 




14.7 






12.279 


juloa 


14.7^* 


1.13 (0.14) 


5.0 


0.95 


(0.12) 




10.5 


GSS 30 


NIRSPEC/Kcck 


4.695 


apr02 


1.2* 


1.97 (0.17) 


8.7 


7.38 


(0.62) 




28.0 


GV Tau N 


TEXES/Gemini 


12.279 


nov06 


25.6'' 


0.55 (0.07) 


3.4 


0.27 


(0.04) 




8.5 


HL Tau 


TEXES/IRTF 


12.279 


dcc02 


7.6° 


1.16 (0.18) 


7.1 


1.82 


(0.29) 




10.9 


HL Tau 


TEXES/Gemini 


8.025 


nov06 


5.5° 


0.34 (0.11) 


2.1 


0.50 


(0.16) 




7.1 




12.279 


nov06 


7.6° 


0.84 (0.13) 


5.1 


1.36 


(0.22) 




12.3 






17.035 


nov06 


13.4° 


0.84 (0.23) 


5.1 


1.07 


(0.29) 




11.6 


HL Tau 


NIRSPEC/Kcck 


4.695 


oct01-nov03 


1.6'' 


0.09 (0.03) 


0.6 


0.27 


(0.06) 




11.9 


In units of 10 ergs s ^ 


cm ^, Value in parentheses is 1 


-a error'' Measured value° ISO SWS archive'' 


Spitzer 



IRS" 3a hmit assuming FWHM = 10.4 km s'^ 



BITNER ET AL. 



TABLE 3 
Results of LTE Model Fits 



Star Telescope One Component Two Component 







Temperature (K) 


Mass (Mg) 


'^cold (K) 


M^old (M®) 


'ihot (K) 


Mhot (M(b) 


AB Aur'' 
GSS 30^ 
HL Tau» 

Elias 29'" 
Elias 29"° 


Gemini 
Gemini 
Gemini 

Gemini 
IRTF 


670 (40) 
535 (45) 
465 (20) 
1210 (90) 
1000 (90) 


0.52 (0.15) 
0.78 (0.16) 
1.09 (0.14) 
0.19 (0.03) 
0.78 (0.07) 


320 (60) 
520 (60) 
460 (20) 


1.65 (0.52) 
0.82 (0.27) 
1.11 (0.17) 


1470 (100) 
3330 (130) 
17QO+'30 

J^""-'_440 


0.076 (0.01) 
002+° °°2 
001+°-°°'' 



^ One component fit to 3(1), S(2), and 8(4). Two component fit includes 8(9).'' One component fit to 8(1), 8(2), 
and S(9). 
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TABLE 4 

Summary of Upper Limits - IRTF/TEXES 



Star 


A 
(/zm) 


Date 




Line Flux'' 


200 K 




Mass in Mjj, 
500 K 


p 




800 K 




51 Oph 


12.279 


jun03 


10.9 




2.2 


< 1.84 






< 3.0 


X 10" 


-2 


< 


1.3 




IQ- 


-2 


AB Aur 


8.025 


oct04 


12.9 




1 9 








<^ 0.0 


V 1 0" 
A XU 


-2 


< 


'\ Q 
o.y 


X 


J-U 


-3 




17.035 


dec02, jan04 


25.2 




1.1 


<^ 4 5 X 


10~ 


-2 


^ ^ 4 


X 10" 


-3 


< 


4.0 




10" 


-3 


AS 209 


12.279 


junOS 


2.5 




0.61 


< 4.2 X 


10~ 


-1 


< 6.8 


X 10- 


-3 


< 


3.1 




10~ 


-3 


AS 353 


17.035 


jun03, jul03 


1.6 




1 .u 


^ V 1 V 




-2 


^ 0.0 


X 10" 


-3 


< 






1 


-3 


FU Ori 


12.279 


jan04, jan05 


6.4 




U.o 


^ Q 7 






^ 1 


X 10" 


-1 


< 


7 1 




lU 


-2 




17.035 


jan05 


5.8 




o.u 


1.0 






^ 1.0 


X 10" 


-1 


< 








-1 


GG Tau 


12.279 


jan04 


1.6 




2.1 


< 2.0 






< 3.2 


X 10" 


-2 


< 


1.5 




IQ- 


-2 




17.035 


novOl 


0.5 




2 8 


< 1 1 X 




-1 


1.0 


X 10" 


-2 


< 


Q 7 






-3 


GW Ori 


12.279 


dccOO, jan05 


8.6 




1 T 
± . 1 


^ lU.O 






^ i. / 




-1 


< 


7 Q 


X 




-2 




17.035 


novOl 


4.9 


< 


Z.o 


< i.i 






^ 1 Q 

< i.o 


X lU 


-1 


< 


1.0 


X 


1 n~ 
iU 


\ 


HD 163296 


12.279 


jun03, jul03 


11.4 


< 


0.6 


< 4.3 X 


10" 


-1 


< 7.0 


X 10- 


-3 


< 


3.2 


X 


IQ- 


-3 




17.035 


jul03 


13.7 


< 


2.35 


< 6.9 X 


10- 


-2 


< 8.3 


X 10- 


-3 


< 


6.2 


X 


10- 


-3 


IRAS 04278+2253 


12.279 


jan04, jan05 


8.0 


< 


0.62 


< 5.9 X 


10- 


-1 


< 9.5 


X 10- 


-3 


< 


4.3 


X 


10- 


-3 




17.035 


jan05 


10.7 


< 


2.8 


< 1.1 X 


10- 


-1 


< 1.3 


X 10- 


-2 


< 


9.7 


X 


10- 


-3 


L1551 IRS 5 


12.279 


dcc02 


13.2 


< 


1.4 


< 1.3 






< 2.1 


X 10- 


-2 


< 


9.7 


X 


10" 


-3 




17.035 


novOO 


20.0 


< 


1.9 


< 7.4 X 


10- 


-2 


< 8.8 


X 10- 


-3 


< 


6.6 


X 


10- 


-3 


Lk Ha 225 


12.279 


jun03 


37.9 


< 


2.5 


< 121.7 






< 2.0 






< 


8.9 


X 


10- 


-1 




17.035 


oct04 


48.2 


< 


3.56 


< 7.1 






< 8.4 


X 10- 


-1 


< 


6.3 


X 


10- 


-1 


SVS 13 


12.279 


dcc02 


13.3 


< 


1.02 


< 4.5 






< 7.2 


X 10- 


-2 


< 


3.3 


X 


10- 


-2 




17.035 


dcc02 


18.8 


< 


2.0 


< 3.6 X 


10- 


-1 


< 4.3 


X 10- 


-2 


< 


3.2 


X 


10- 


-2 


V892 Tau 


12.279 


jan04 


31.2 


< 


3.0 


< 2.86 






< 4.6 


X 10" 


-2 


< 


2.1 


X 


10- 


-2 




17.035 


dec02 


71.7 


< 


4.9 


< 1.9 X 


10- 


-1 


< 2.2 


X 10" 


-2 


< 


1.7 


X 


10- 


-2 


V1057 Cyg 


12.279 


oct04 


7.0 


< 


0.8 


< 14.0 






< 2.3 


X 10" 


-1 


< 


1.0 


X 


10- 


-1 


Z CMa 


12.279 


dec03 


142.6 


< 


3.7 


< 238.3 






< 3.8 






< 


1.74 








17.035 


dec02 


170.8 


< 


14.0 


< 36.8 






< 4.4 






< 


3.27 







* Measured TEXES flux value'' 3(t limit in units of lO"!" ergs s"! cm-2 
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TABLE 5 

Summary of Upper, Limits - Gemini/TEXES 



Star 


A 


Date 








Mass in Mj„p 








(/im) 




(Jy) 


Line Flux'''*' 


200 K 


500 K 


800 K 


49 Ceti 


12.279 


nov06 


0.2^^ 


< 0.17 


< 3.1 X 10 


< 5.0 X lU 


< 2.2 X 


10—* 




17.035 


nov06 


0.19"^ 


< 0.88 


< 6.5 X 10 


< 7.7 X 10 


< 5.8 X 


10—* 


AS 209 


12.279 


jul06 


2.0"* 


< 0.34 


< 2.3 X 10"* 


< 3.8 X IQ-'^ 


< 1.7 X 


10-3 




17.035 


jul06 


4.4° 


< 1.03 


< 2.9 X 10"^ 


< 3.4 X IQ-'^ 


< 2.6 X 


10-3 


AS 353 


12.279 


jul06 


1.11' 


< 0.29 


< 3.2 X 10~* 


< 5.1 X IQ-'^ 


< 2.3 X 


10-3 


DoAr 21 


12.279 


jul06 


0.14^ 


< 0.18 


< 8.4 X 10~* 


< 1.4 X 10-2 


< 6.1 X 


10-3 


Elias 29 


17.035 


jul06 


30.5" 


< 1.36 


< 3.8 X 10~^ 


< 4.6 X IQ-'^ 


< 3.4 X 


10-3 


FU Ori 


8.025 


nov06 


3.5^ 


< 1.14 


< 5143.6 


^ O O ^ > "in 1 

< 3.8 X 10 


< 4.5 X 


10-2 




12.279 


nov06 


3.6^ 


< 0.22 


< 2.7 


< 4.3 X 10 


< 2.0 X 


10-2 


GM Aur 


12.279 


nov06 


0.25'' 


< 0.54 


< 5.2 X 10 


< 8.3 X 10 ^ 


< 3.8 X 


10-3 


Gbb 30 


17.035 


jul06 


30.75* 


< 1.3 


< 3.7 X 10 


< 4.3 X 10 ^ 


< 3.3 X 


10-3 


HD 1415d9 


12.279 


jul06 


1.13" 


< 2.75 


< 1.3 


^ o "1 ^ , 1 n — 2 

< 2.1 X 10 


< 9.6 X 


10-3 


tiD loo2yo 


12.279 


jul06 


14.2" 


< 1.18 


< 8.6 X 10 


. -1 A in — 2 

< 1.4 X lU 


< 6.2 X 


10-3 




17.035 


jul06 


21.6" 


< 2.42 


< 7.2 X 10~^ 


^ s V 1 

\ O.O A _LU 


< 6.4 X 


10-3 


LkCa 15 


12.279 


nov06 


0.12S 


< 0.08 


< 7.6 X 10-2 


< 1.2 X 10-3 


< 5.6 X 


10-* 




17.035 


nov06 


0.48" 


< 0.97 


< 3.8 X 10-2 


< 4.5 X 10-3 


< 3.3 X 


10-3 


MWC 758 


17.035 


nov06 


3.8" 


< 0.65 


< 5.2 X 10-2 


< 6.1 X 10-3 


< 4.6 X 


10-3 


RW Aur 


8.025 


nov06 


0.07* 


< 0.37 


< 130.9 


< 9.7 X 10-3 


< 1.1 X 


10-3 




17.035 


nov06 


1.87* 


< 0.29 


< 1.1 X 10-2 


< 3.7 X 10-2 


< 1.0 X 


10-3 


TW Hya 


12.279 


feb06 


0.5* 


< 0.6 


< 7.6 X 10-2 


< 1.2 X 10-3 


< 5.5 X 


10-* 




17.035 


feb06 


1.4' 


< 0.7 


< 8.9 X 10-2 


< 1.4 X 10-3 


< 6.5 X 10-"* 


V 1331 Cyg 


17.035 


nov06 


1.65* 


< 0.52 


< 3.1 X 10-* 


< 3.7 X 10-2 


< 2.8 : 


X 10-3 


VV Ser 


12.279 


jul06 


4.61*' 


< 0.55 


< 1.6 


< 2.6 X 10-2 


< 1.2 X 10-2 




17.035 


jul06 


2.6f5 


< 1.2 


< 1.4 X 10-* 


< 1.7 X 10-2 


< 1.3 X 10-2 



^ 3<T limit in units of 10 ** ergs s * cm 2b Upper limits calculated assuming FWHM=5.5 km s * for TW 
Hya, 10 km s-* for others" rWahhai et al.l l l2007f )'* ILiu et al.l llT996t) " ISO SWS archive* Spitzer IRSS Measured 
TEXES fiux value*' IRAS* IRatzka et aTl | |2007D 
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Fig. 13. — Upper limits for IRTF H2 S(l) observations near 17 //m. The overplotted Gaussian is centered at the stellar velocity of each 
source and represents the 3-a upper limit based on an assumed FWHM of 10 km s~^. 
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Fig. 14. — Upper limits for IRTF H2 S(2) observations near 12 //m. The overplotted Gaussian is centered at the stellar velocity of each 
source and represents the 3-a upper limit based on an assumed FWHM of 10 km s~^. 
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Fig. 15. — Upper limits for Gemini H2 S(l) observations near 17 /xm. The overplotted Gaussian is centered at the stellar velocity of each 
source and represents the 3-a upper limit based on an assumed FWHM of 10 km s~^. 
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Fig. 16. — Upper limits for Gemini H2 S(2) observations near 12 /im. The overplotted Gaussian is centered at the stellar velocity of each 
source and represents the 3-a upper limit based on an assumed FWHM of 10 km s~^. 
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APPENDIX 
ACCRETION SHOCK HEATING 
The total luminosity L from a shock of area A is given by: 

L = (Al) 

where Po is the preshock mass density and Vg is the shock velocity (the normal component of the velocity of the flow 
with respect to the shock surface). The mass accretion rate Mace through the shock is given by: 

Mace = PoVgA = mpHoVsA, (A2) 

where Hq is the gas hydrogen nucleus number density and is the mass per hydrogen nucleus 2.3 x 10^^"* gm)- 
Therefore, the luminosity can be rewritten 

L = ^Maccvl (A3) 

However, in an accretion shock onto an optically thick disk, where one half of the radiation is em itted toward the disk 
midplane and is absorbed, the escaping luminosity is given by Csee lNeufeld fc HollenbachlfT99l : 

L^sk = ^Maccvl (A4) 

The shock velocity is on the order of (but somewhat smaller due to the oblique incident angle of the infall to the shock 
front) the freefall velocity onto the disk, or, for our r > 10 AU constraint, ~ 5 — 10 km s~^. Assuming the accretion 
rate from the core onto the disk through the accretion shock is similar to the accretion rate from the disk onto the 
star, the measured accretion rates for our sources are of order Mace ^ 10^^ Mq yr^^. Therefore, 

Ld^,k^ 4.x 10-^ M^rv% Lq, (A5) 

where M_7 = M© yr^^ and Vsq = cm s"i== Vg/lO km s~^ 

To determine the luminosity in the pure rotational lines of H2 , one needs the fraction fj of the total shock luminosity 
that emerges in the H2 0-0 S(J) line. This fraction depends on the preshock density rip, the sho ck velocity Vs, the 
am ount of depletion of the pre shock gas coolants, and whether the shock is "C type" (iDraind 1980) or "J type" 
(cf. iHollenbach fc McKedHOTOl ). The preshock density can be estimated by taking the shock area A to be at least 
27rrg ~ 1.4 X 10^^ cm^, with ~ 10 AU and the factor of 2 to account for both sides of the disk. Using Eq. (2), we 
obtain 

no ~ 3 X lO^A/.TW^giAgo^ cm^^, (A6) 

with A30 = A/lO^o cm2. iBurton et all (|1992l ) present results for J shocks with Uo — 10 cm ^ and Ws = 5 — 10 km s 
in terms of the intensity /,/ of an H2 0-0 S(J) line. Here, 

= TT^ —■ 

u.ompnoVs 

They show cases with high abundances of gas phase oxygen not in CO and with extremely low abundances of gas 
phase oxygen not in CO. The latter case is perhaps more realistic since oxygen not in CO is expected to freeze out 
as water ice in dense cores , and a slow 5-10 km s~^ shock does not release water from the ice mantles to the gas 
(|Hollenbach fc McKeelll979f ). In the former case, the gas phase oxygen not in CO rapidly co nverts to gas phase H 2O 
in the shock, and the H2O dominates the shock cooling and thereby weakens the H2 lines. iBurton et al.l (|1992l see 
Figures 5b and 6) find that for the case with H2O freezeout /i ^ 2.5 x 10~^, /2 ~ 2.5 x 10~^ and fs ^ 10~^ . In the 
case with abundant gas phase oxygen not in CO, /i ~ 10"^, /2 ~ lO"'*, and /a ^ 10~'^. IBurton et al.l (jl992l ) do not 
present results for f^, but inspection of the unpublished output of these runs reveals ~ O.5/3. A lower results in 
lo wer ratios oi JaI f2- 

iDraine et al.l ([1983) present /,/ for C type shocks with Uq = 10^ cm^^, preshock magnetic field component parallel to 
the shock front Bo = 0.5 mG, electron abundance Xc — 10^^, and where gas phase oxygen not in CO is not depleted. 
Their result at Vg — 10 km s^^ implies /i ~ 10~^, /2 4 x 10^^, and /a ^ 10~^. They also do not present results for 
/4, but we again estimate /4 O.5/3. Presumably, these fractions would be somewhat larger if oxygen is allowed to 
freezeout as water ice. 

From the above results, we see for J shocks with water ice freezeout and for C shocks, the luminosities in the H— 2 
0-0 S(l), S(2) and S(4) lines are ~ 10~® — 10~^M_7W^g Lq. These luminosities (assuming Af_7 ~ 1 correspond to 
our observed luminosities, suggesting accretion shocks as a viable excitation mechanism for these emission lines. For 
J shocks without freezeout the luminosities are about an order of magnitude less than observed for Af_7 ^ 1. 
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The se same shock resuhs can be used to estimate the intensities of CO mid J transitions in the shocks. iBurton et al.l 
(|1992| ) show that in the J sh ocks p resented here the CO J= 17-16 transition is '--^ 3 times stronger than the H2 0-0 
S(l) hne. [D raine fc Robergd (|1984r ) and lDraine et alj (|1983f ) present C shock models from which the CO J— > J-1 hne 
strengths can be estimated. For Uq = 10^ cm~^, Vg = 10 km s~^, and Bq and from above, the CO J J-1 hne 
intensities peak at J ~ 10 — 12 with strengths comparable to that of H2 0-0 S(l). Thus, Herschel observations of these 
disks with the HIFI instrument may detect the submiUimeter CO hnes near the J peak, vaUdate the shock models, 
and constrain the shock parameters such as rio, Vs, and A. 

It is difficult to estimate the strengths of lines other than those of H2 and CO because their strength depends on 
their gas phase abundances, and all species other than the undepleted H2 and CO (e.g., Fe, Fe"^, S, O) may be heavily 
depleted on grains as ice mantles and/or refractory material. Such slow shocks are unlikely to remove them. Radiative 
transfer calculations need to be performed to see if the dust grains are warm enough to thermally sublimate the ice 
mantles at distances of 10-30 AU from the star. A more detailed shock modeling of some of these observed sources is 
warranted. 



